Athermal Shear-Transformation-Zone Theory of Amorphous Plastic Deformation I: 

Basic Principles 
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We develop an athermal version of the shear-transformation-zone (STZ) theory of amorphous 
plasticity in materials where thermal activation of irreversible molecular rearrangements is negligible 
or nonexistent. In many respects, this theory has broader applicability and yet is simpler than its 
thermal predecessors. For example, it needs no special effort to assure consistency with the laws 
of thermodynamics, and the interpretation of yielding as an exchange of dynamic stability between 
jammed and flowing states is clearer than before. The athermal theory presented here incorporates 
an explicit distribution of STZ transition thresholds. Although this theory contains no conventional 
thermal fluctuations, the concept of an effective temperature is essential for understanding how the 
STZ density is related to the state of disorder of the system. 



I. INTRODUCTION 

The shear-transformation-zone (STZ) theory of amor- 
phous plasticity, to date, has been applied most success- 
fully to "thermal" glassy systems at temperatures high 
enough that they exhibit linear viscosity and that non- 
linear flow at larger driving stresses is controlled by ther- 
mally activated processes 0, Q- Our purpose here is 
to examine the opposite, "athermal" situation, where 
the ambient temperature is negligible, and all rearrange- 
ments of the constituent elements are driven entirely by 
applied forces. Systems of the kind to be discussed here 
include noncrystalline solids well below their glass tem- 
peratures, dense granular materials, and various kinds of 
soft materials such as foams, colloids, and the like. 

An increasingly useful source of information about 
these systems is numerical simulation which, while lim- 
ited in comparison with laboratory experiments on real 
materials, has certain compensating advantages. For ex- 
ample, athermal materials seem to be intrinsically un- 
stable against nonuniform failure via shear banding or 
fracture. Such instabilities are much more difficult to 
observe and control in the laboratory than in large-scale 
computations. One of our main goals in this project is 
to develop a predictive description of athermal plasticity 
- an analog of the Navier-Stokes equation for amorphous 
solids - that can be used in heterogeneous situations. In 
the present paper and its sequel |3|, however, we confine 
our attention to spatially homogeneous systems and test 
our results by comparing with simulations rather than ex- 
periments. Another advantage of numerical simulation is 
that it allows us to observe internal states of the system 
that are not easily accessible in laboratory experiments. 
That capability is a central feature of the following paper. 

From its inception, the STZ theory has been built 
on the flow-defect theories of TurnbulL Cohen, Argon, 
Spaepen, and others. USEIlllllllEEllt describes 
plastic deformation in amorphous solids, or solidlike ma- 
terials, but not in liquids. The assumption is that ir- 
reversible molecular (or granular) rearrangements occur 



only at sparsely distributed sites - the STZ's - within 
an otherwise elastic material. The validity of this as- 
sumption was demonstrated explicitly in | l2| , but it goes 
back to essentially all of the previously cited earlier work. 
The STZ model is strictly valid only when the local rear- 
rangements occur infrequently and independently of each 
other, and when they require either substantial thermal 
activation or, in the athermal situations of interest here, 
sufficiently large external driving forces. If the activa- 
tion energy or work needed to drive a rearrangement is 
small of order ksT, or if the sites at which rearrange- 
ments occur cover most of the system, then the material 
is effectively a liquid and STZ theory is not applicable. 

We visualize an STZ as a localized group of molecules 
that is more susceptible than its neighbors to a shear- 
ing transformation in some direction. That is, these 
molecules must collectively surmount only a relatively 
small energy barrier in order to undergo an irreversible 
shear. Once this happens, it seems reasonable to sup- 
pose that they will resist further shear in the original 
direction, but may be especially susceptible to a reverse 
shear. One might imagine that the first transition has 
redistributed the local stresses in such a way as to fa- 
vor a reverse transition if the applied stress changes sign. 
We see no strong requirement that the reverse transition 
must bring the molecules back to exactly their original 
positions; but it is this approximate picture that suggests 
a two-state model of STZ's. 

A primary rationale for the two-state model is that it 
provides a simple mechanism by which the system re- 
tains orientational memory of prior deformations. Along 
with the dynamic transition between jammed and flow- 
ing states, orientational memory is one of the universal 
features of amorphous plasticity that we believe must be 
captured by any satisfactory theory. A related require- 
ment for an acceptable theory is that it must include a 
mechanism by which orientational memory is lost during 
deformation. Here that mechanism is the annihilation 
and creation of STZ's at a rate proportional to the rate 
of energy dissipation. Under athermal conditions, anni- 
hilation and creation occur only in response to STZ tran- 
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sitions; thus this mechanism may also be seen as a rough 
description of the cascades of rearrangements following 
STZ-like events seen by Maloney and Lemaitre 0, , 
and emphasized by Argon and Demkowicz in papers to 
be discussed in a sequel to this one. [jjj, Uji Il3 

The defining feature of an athermal system is the con- 
straint that, because thermal activation of transitions is 
negligible or nonexistent, molecular rearrangements oc- 
cur only in response to driving forces. No motion oc- 
curs in the absence of such forces, and no rearrangement 
moves in a direction opposite to that in which the force is 
applied. In other words, molecular configurations cannot 
move uphill in energy as they may when thermal fluctu- 
ations are present. Stress-induced shear transformations 
are intrinsically irreversible events. Work is done on the 
system as the STZ's are driven over energy barriers, and 
energy is dissipated as the system moves downhill toward 
new stable states. In contrast to the picture proposed in 
earlier papers 0, |2j| |2(| which assumed only one kind 
of STZ, the model to be developed here allows STZ's to 
occur with a range of different sizes and transition thresh- 
olds. With this generalization, the STZ model may un- 
dergo limited irreversible deformations when the applied 
stresses are less than the nominal yield stress. The onset 
of athermal flow at an apparent yield stress then must be 
a dynamic phenomenon. As in the original STZ theories, 
it occurs when there is an exchange of stability between 
the jammed steady states, where nothing is moving, and 
the flowing steady states, where motion-induced annihi- 
lation and creation of STZ's balances the rate at which 
zones are inactivated by forward transitions. In this way, 
the STZ picture describes the dynamics of plastic yield- 
ing by the same mechanism that it uses to describe the 
memory effects mentioned above. 

Although ordinary thermal fluctuations are absent in 
the athermal models discussed here, the concept of an 
effective disorder temperature is essential. Some of the 
earliest work in this field recognized that the density of 
flow defects could be related to an intensive quantity 
such as the free volume vj (the inverse of the deriva- 
tive of an entropy with respect to the volume). Inten- 
sive quantities of this kind characterize the state of the 
system as a whole and not just that of a subset of its 
degrees of freedom. Thus Cohen and Turnbull in 1959 
(see also Spaepen p|) proposed that the density of 
flow defects in an amorphous solid be proportional to 
cxp (— const. /vf), and not just to Vf itself. In 2|, one of 
us (JSL) argued that the appropriate generalization of 
free volume is an effective temperature T e ff that char- 
acterizes the state of configurational disorder in the sys- 
tem. T e ff equilibrates to the ambient temperature T 
at high T, but may fall out of equilibrium at low T 
where disorder is generated by the atomic-scale, configu- 
rational rearrangements that accompany mechanical de- 
formation. If T e ff — (EsTz/ks) X, where E$tz is a 
characteristic STZ formation energy, then the STZ den- 
sity is proportional to the Boltzmann factor exp (— 1/x)- 
This is a direct analog of the free-volume formula and, 



in fact, reduces to it in the case of a system under con- 
stant pressure with a positive "effective" thermal expan- 
sion coefficient. Importantly, the time variation of the 
STZ density is governed by the dynamics of x- (F° r 
more information about effective-temperature theories, 

see [n us m m h m . ) 

The scheme of this paper is as follows. We review and 
reformulate basic features of the STZ theory in Section 
II. Then, in Section III, we derive the athermal equations 
of motion for the STZ state variables, propose a specific 
form for the athermal rate factor, and show a few illus- 
trative examples of how the theory behaves in various 
experimental situations and with various choices of the 
constituent parameters. 

II. STZ BASICS 

The STZ theory is a phenomenological construction. 
Our strategy has been to start with what amounts to a 
caricature of an amorphous material, specifically, a model 
in which applied stresses and two-state STZ's remain 
aligned along fixed axes. We deduce from this rudimen- 
tary model the internal state variables that are needed 
to describe its behavior, and then derive equations of 
motion for those variables. When applying the theory 
to more realistic situations in three dimensions, at least 
in simple geometries such as those we encounter here, 
we assume that we can retain the form of our equations 
of motion but replace certain state variables by tensors 
when required by symmetry. In short, we see how far we 
can go with minimal models and, to the extent possible, 
test these models by comparing our theoretical predic- 
tions with experimental data as in Q, Q • 

Accordingly, we start by considering a two-dimensional 
system, and subject it only to pure shear deformation ori- 
ented along a fixed pair of principal axes, x and y. It is 
sufficient for present purposes to assume that the popula- 
tion of STZ's consists simply of zones oriented along the 
two principal axes of the deviatoric stress tensor, which 
we take to be s xx — —s yy — s and s xy = 0. Choose the 
"+" zones to be oriented (elongated) along the x axis, 
and the "— " zones along the y axis. We assume that the 
STZ's occur in many different varieties, with the sym- 
bol a representing, for example, their actual orientations 
with respect to the stress axes as in 26] or, explicitly in 
what follows, their transition thresholds. Thus we denote 
the population density of zones oriented in the "+/— " di- 
rections by the symbol n±(cn). 

With these conventions, the plastic strain rate - more 
generally, the plastic part of the rate-of-deformation ten- 
sor - is: 

D pl = —D pl = D pl 

xx yy — 

= — J da(R a {~s) n-(a) - R a {-s) n+(a)) . (2.1) 

Here, A is a material-specific parameter with the dimen- 
sions of volume (or area in strictly two-dimensional mod- 
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els), which must have roughly the same order of magni- 
tude as the volume of an STZ, that is, a few cubic or 
square atomic spacings. tq sets a time scale for these 
processes. The integration is over the relevant space of 
parameters a. The integrand is the net rate per unit 
volume at which a-type STZ's transform from "— " to 
"+" orientations. R a (s)/To and i? a (— s)/to are the rates 
for forward ("-" to "+") and backward ("+" to "-") 
transitions respectively. For later convenience, we have 
written these rates as functions of a dimensionless stress 
5 = s/sy, where s y will turn out to be the dynamic yield 
stress in the athermal theory. At present, we need to 
think of s y only as a characteristic scale for measuring 
stresses. 

The next step is to postulate a master equation for the 
populations n±(a). As before, we do this in a mean- field 
approximation. Using earlier notation, we write 

Toh±(a) — R a (±s) n T (a) — R a (=F§) n±(a) 

+ m (^Me-^-n ± (a^j. (2.2) 

The first pair of terms on the right-hand side describes 
the same switching back and forth of the STZ's that ap- 
pears in Eq. (|2.1|) .The last terms describe the creation 
and annihilation of zones at a rate T(s)/tq that, in a 
mean-field sense, we assume to be the same for both 
n+(a) and n_(a), independent of a and of local prop- 
erties of the system. T is a non-negative, scalar quan- 
tity that vanishes when the rate of deformation is zero; 
in earlier papers |25l |26( we have argued that it must 
be proportional to the rate per STZ at which mechan- 
ical work is dissipated via irreversible plastic deforma- 
tion, exp (— 1/x) is the steady-state density of STZ's 
achieved by the system during persistent deformation. 
As discussed in the Introduction, x 1S the effective tem- 
perature measured in units of the STZ formation energy. 
In writing this part of Eq. I|2.2|l . we are using the princi- 
ple of detailed balance to fix the ratio of the annihilation 
and creation rates. Note that Eq. (|2.2|) contains no aging 
or spontaneous relaxation, consistent with the assump- 
tion that thermal fluctuations are absent. 

Now suppose that n±(a) = n±p(a) and n x (a) = 
noop(a), where p{a) is a normalized distribution over 
a. At this point, we make the important simplifying as- 
sumption that p(a) is not itself a dynamical quantity that 
changes during deformation. Performing the integration 
in Eq. l|2.1[l . we find that 



where 



D pl = - (r(s) n_ - R(-s) ?i+ 
to V 



-^(s) = / R a (s) p(a) da. 



(2.3) 



(2.4) 



Similarly, we integrate both sides of Eq. (|2.2II over a to 
obtain 



tq n± = R(±s) n T — R{Ts) n± + T(s) (- 



■ e — n± 



Thus we recover exactly the earlier equations of motion 
for the STZ populations, but with a modified interpreta- 
tion of the rate factors. 

Before writing an equation of motion for \, and then 
moving on to the specifics of the athermal theory, we rein- 
troduce some convenient notation and rewrite the equa- 
tions of motion in the form in which we shall use them. 
We define the following dimensionless quantities: 



6q = Arioo , 

n + + 7i- 
A = — 



n + — ri_ 
n + + n _ 



C(S) ee - (R(s) + R(-S)) , 



T(S) 



R(s) - R(-s) 



(2.6) 
(2.7) 

(2.8) 

(2.9) 

(2.10) 



R(s) + R(-s) ' 
Then Eq. (£21 is 

T DP l = e AC(s)(T(s)-m); (2.11) 
and Eq. I|2.5I) becomes a pair of equations for m and A: 



T 777 = 2 C(S) (T(af) - m) - e~^> 



and 



r A - r(S) (e 



-i/x 



A). 



(2.12) 



(2.13) 



Using the preceding notation, we choose the equation 
of motion for the effective temperature \ t° be the ather- 
mal version of Eq. (3.5) in [2J. So far, none of the in- 
gredients of our equations of motion have been written 
specifically in their athermal forms; but here we deviate 
by dropping terms that refer explicitly to mechanisms by 
which x relaxes to the ambient temperature. Thus: 



tqCqx = eo AT(s) (xoo - X)- 



(2.14) 



(2.5) 



This is basically an expression of the first law of thermo- 
dynamics. The left-hand side is the time rate of change of 
the configurational internal energy roughly approximated 
as the product of the effective temperature x multiplied 
by a specific heat cq. The latter quantity is expressed 
in units of fee per atom and thus is of order unity. The 
right-hand side of Eq. Q2.14|) is proportional to the rate 
of energy dissipation per unit volume, that is, the dissi- 
pation rate per STZ, T multiplied by the STZ density, 
A. 

The last factor on the right-hand side of Eq. (|2.14l) ap- 
pears because there must be an upper bound on x'i the 
disorder temperature cannot simply increase indefinitely 
under continued athermal deformation but, rather, must 
settle to some steady-state value Xoo- That limiting be- 
havior has been demonstrated explicitly by Ono et al. in 
simulations of a sheared foam |l9j , in which the authors 
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showed that a variety of different definitions of an effec- 
tive temperature for an athermal system are consistent 
with each other. They also found that Xoo has a nonzero 
value at sufficiently small strain rates and increases, as 
argued in [2J, only when the strain rate becomes compa- 
rable to other relevant, internal, inverse time scales. In 
0, the small-strain-rate value of Xoo was estimated to 
be roughly the ratio of the yield stress to the shear mod- 
ulus or, more or less equivalently, the ratio of the glass 
temperature to the STZ formation energy. If the former 
estimate can be taken literally, then Johnson's analysis 
of yield strengths in a wide range of glasses in [271 im- 
plies that Xoo is a universal number of order 0.02-0.04. 
That value is consistent with the one found in Q , where 
a direct estimate of the STZ energy in a metallic glass 
was available from viscosity measurements. 

The preceding estimates of Xoo give us useful insight 
regarding the general structure of the STZ theory sum- 
marized by Eqs. H2.11J1 through H2.14[l . In our atomic 
units, the density should be of order unity, and eo 
as defined in Eq. H2.6(l also must be of order unity. 
Thus, if Xoo is small of order 0.1 or less, the density 
of STZ's, A = exp(-l/xoo), is of order 10~ 3 or appre- 
ciably smaller, consistent with our basic assumption that 
the STZ's are rare defects that interact only weakly with 
each other. In retrospect, we recognize that earlier STZ 
theories that did not include the effective temperature, 
e.g. [ij, required improbably large values of the STZ 
density in order to agree with experiment. 

A appears as a rate-determining prefactor on the right- 
hand sides of Eqs. I|2.11|l and (|2.14|) . which govern the 
bulk system-wide variables D pl and \\ but A does not 
appear in a similar way in Eqs. H2.12[l or (|2.13|) . which 
pertain to the dynamics of individual STZ's. It follows 
that the plastic strain rate and the effective temperature 
respond much more slowly to changes in stress than the 
do the internal STZ variables m and A, and that the 
slow dynamics of the effective temperature controls the 
observable mechanical behavior of the system in most 
circumstances. 



III. ATHERMAL THEORY 

The crux of the athermal STZ theory is the choice 
of the rate factor R(s) and the resulting expression for 
the creation and annihilation factor T. An immediate 
and important simplification follows from the athermal 
constraint that no motion occurs in a direction opposite 
to that of the applied force. Thus R(s) must vanish if 
s < 0. We then find from Eq. (|2~TU|) that 

T(s) = sign(s) = 1 I T . (3.1) 
\s\ 

This result is trivially correct no matter how complicated 
the transition rate might be. It immediately tells us that 
jamming - i.e. D pl — - occurs only for m — ±1, de- 
pending on the sign of the stress, and that no transitions 



in either the forward or backward direction are occur- 
ring in the jammed state. For example, when s > in 
Eq. I|2.3[) . jamming occurs because both n_ = and 
R(-s) = 0. 

A second immediate simplification comes from the ob- 
servation that, with no uphill transitions in energy, all 
the work done on the system must be dissipated and 
none can be stored internally. Thus the dissipation rate 
per STZ, in the dimensionless form required by Ea. l|2.5|l . 
must be 

r(a,m)= 2T ° s f Pl = 2sC(s) f^r-mV (3.2) 
eo A \\s\ / 

We thus recover the expression for T originally guessed in 
|12| but now without the sign problem pointed out there. 
The work done by the external driving force cannot be 
negative in the athermal limit because the plastic flow 
must have the same sign as the stress. rEa. (|3.2|l can 
be derived systematically using Pechenik's method [2^, 
[2(J. That analysis reveals that, in the athermal limit, 
the recoverable internal energy vanishes for m 2 < 1.) 

To complete the definition of our model, we must spec- 
ify a form for R(s). Consider just a single species of 
STZ, and suppose that the parameter a determines the 
activation threshold, say, s a . That is, let R a (s) vanish 
for s < s Q , with the understanding that the absence of 
thermal fluctuations means that the only allowed tran- 
sitions are those that are driven mechanically by suffi- 
ciently large stresses. A convenient, minimal form for 
R a (s) is 

R a (s) = /2(s-Sa) fors>s a (33) 
w [0 for -oo < s < s a . y ' 

This expression is a rough description of an athermal sys- 
tem moving in a double-well potential. The system re- 
mains trapped until the applied force s reaches its thresh- 
old Sq, , at which point the unbalanced force rises linearly 
in s — s Q . If the response is dissipative, i.e. frictional 
or viscous, then the speed at which the system moves 
away from its original position will also be proportional 
initially to s — s a . In our case, the proportionality coef- 
ficient is incorporated into the factor Tq 1 . 

For purely athermal systems such as granular materials 
or foams where there are no other relevant time scales, 
this interpretation of the rate factor makes sense only if 
To is comparable to, or longer than the inverse of the total 
strain rate. In the opposite limit, where the duration of 
an STZ transition is very short compared to the interval 
between transitions, the only relevant time is the inverse 
strain rate itself. Therefore, the rate of STZ transitions 
Tq 1 must be proportional to the total strain rate. This 
is the common limiting situation in which the number of 
irreversible atomic rearrangements, i.e. STZ transitions, 
does not depend on the length of time during which the 
material has been loaded but only on the extent of the 
deformation. The situation is quite different in molecu- 
lar solids, even at athermally low temperatures, because 
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then the molecular vibration frequency governs the rate 
of the forward, stress-enabled transitions across the STZ 
energy barrier. 

According to Ea. (|2.4|l . R(s) is the average over rate 
factors R a (s) with a normalized weight factor that we 
can denote by p(s a ). In choosing p(s a ), we are led by 
the following considerations: 

• If characteristic stresses for a system of interest are 
of order our scale factor s y , and if s y is the only 
stress scale in the problem, then p(s a ) should have 
a peak at s a = 1. 

• The probability of very low thresholds, s„ < 1, 
must be vanishingly small. The athermal con- 
straint means that there are no negative thresh- 
olds, i.e. no transitions in the direction opposite to 
the stress. Therefore p(s a ) must be such that R(s) 
vanishes smoothly at s = and remains zero for all 
s < 0. 

• We expect that different materials, under different 
circumstances, will have different threshold distri- 
butions, and we therefore need at least one extra 
parameter that controls the width and/or shape 
of p(s a ). Our minimal, phenomenological rule for 
model building implies that we should start by in- 
troducing only one such parameter, denoted below 

byC 

A distribution that satisfies all these criteria and is con- 
venient for numerical purposes, is: 



C c+1 - c 




§i exp(-Cs a ) 



(3.4) 



This distribution has a width of order £ 1 / 2 near its peak 
at s a = 1. R(s) is now obtained from Eq. I|2.4[l : 



c c+1 r 

R(s)=2—^- (s - s a ) exp(-£ s a ) ds a (3.5) 
lo 







from which we see that 
R(s) . 



s c+2 for s -> +0; 
2 (s - 1) for s > 1. 



(3.6) 



Our equations of motion are now conveniently written 
in the form: 



where 



D pl (s,m,A) = —Aq(s,m), 

TO 



q{s,m) = C(s) (j^j- - to); 



m 



l(s,m) (l 



-i/x 



A 



(3.7) 



(3.8) 



(3.9) 



A = — s q(s, to) (e 

TO 



-i/x 



A). 



and 



X 



2e 

CO T 



Asq(s,m) (xoc - x)- 



(3.10) 



(3.11) 



Eqs. I|3.10|l and (|3.11|) tell us that the only stable, 
steady-state solutions of the preceding equations must 
have x — Xoo and A = exp(— 1/xoo), consistent with our 
expectation that the system must flow to a fixed point 
with an STZ density exp(— 1/xoo)- Then Eq. I|3.9|l 
has two stationary solutions: the jammed state with to = 
±1 (depending on the sign of s) and D pl = 0; and the 
flowing state with to = 1/s and D pL ^ 0. These two 
solutions coincide at to = s = ±1. It is easy to check, 
as in many earlier STZ papers, that the jammed state 
is dynamically stable for \s\ < 1 and the flowing state 
for \s\ > 1. An exchange of stability occurs at s y = 1, 
which justifies our earlier choice of units for the stress. 
Note that these conclusions are entirely independent of 
the rate factor, whose stress dependence enters only via 
the function C(s) in the equations of motion. 

To look at the time-dependent behavior of these equa- 
tions, we must include elastic as well as plastic responses. 
Assume that the total rate of deformation D tot is a linear 
superposition of elastic and plastic parts, that is 



D tot = S +D pl(~ S)mjA)i 

2 n 



(3.12) 



where fl is the shear modulus measured in units s y . In 
one common class of experiments, the material is sheared 
at a fixed rate D tot =7/2, and the stress is measured as 
a function of the strain 7. To model such experiments, 
we write Ea. (|3.12|) in the form 



ds / 2e A 

T = V 1 : 9( s > m ) 



(3.13) 



and similarly transform Eqs. I|3.9|) . I|3.10|l and (|3.11l) : 

(3.14) 



dm 2 / tos _ x , 

— = - — q(s,m)[l — e /x 

d'y 7 tq V A 



dA 



— = - — sq(s,m) (e 
«7 7 T o 



-i/x _ 



a); 



and 



dx = 2e 
d'y c 7 r 



Asq{s,m) (xoo - x)- 



(3.15) 



(3.16) 



To illustrate the predictions of this theory, we show 
in Fig. a sequence of stress-strain curves, 5(7), for 
different values of 7 To. A corresponding set of graphs 
of x{l) is shown in the lower panel of that figure. For 
each of these curves we have used £ = 1 (implying a 
broad distribution of low- lying thresholds), \oc = 

1 (a 
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FIG. 1: Upper panel: The normalized stress s/s y as a func- 
tion of strain 7 for different values of the normalized strain 
rate 2D tot ro = 7 to. The parameters used in integrating Eqs. 
ll3~T31l - ll3~ToTl are: £ = 1, Xoo = 1, Xo = 0.5, p, = 45, e = 1 and 
co = 0.25. The different values of the normalized strain rate 
(from top to bottom) are: 2D tot r = jr = 0.3,0.1,0.015. 
The initial values used are: 771(7 = 0) = (i.e. no previous 
deformation) and A(7 = 0) = exp(— 1/xo)- Lower panel: The 
corresponding curves of \ as a function of 7 for the three 
cases. See text for discussion. 



large value, chosen here for illustrative purposes), and 
an initial value of x — Xo = 0.5. We also choose Jx = 45, 
60 = 1) and Co = 0.25. Note that, for small 7 To = 
0.015 (lower curve, upper panel), the flow stress s/ at 
large 7 is approximately equal to the yield stress, s/ = 1, 
and that there is substantial yielding at smaller stresses 
because £ is small. For larger values of 7 To , the stress 
rises nearly elastically to a peak as \ an( i the number 
of STZ's increases, and then drops as the plastic flow 
induces strain softening. The flow stresses are higher in 
these situations. In all cases, the effective temperature \ 
ultimately reaches its steady-state value Xoo- 

A second common class of experiments is that in which 




y 

FIG. 2: The normalized stress s/s y as a function of the strain 
7 for the stress controlled loading shown in the inset. The 
parameters used are the same as those in Fig. Note 
that, although the maximum absolute normalized stress is 
s/s^l = 0.9, significant sub-yield plastic deformation is visi- 
ble, in addition to memory effects. 



the stress rather than the strain is controlled. In this 
situation, we must solve Eq. (|3.12|) as written, with 
jjtot _ on t ue left-hand side and s a predetermined 
function of time t on the right. Figure |2 illustrates a 
stress-strain curve for a case in which the stress is cy- 
cled as shown in the inset. We have chosen the stress 
to remain always less than the yield stress in order to 
illudtrate the effects of small ( = 1. The material pa- 
rameters are the same as in Fig. ^ Note the appearance 
of sub-yield deformation as before, and also the hysteresis 
associated with energy dissipation during that deforma- 
tion. 



IV. CONCLUDING REMARKS 

The athermal STZ theory appears to be cleaner and 
more broadly applicable than its predecessors. It does 
not, of course, replace the thermal STZ theory that is 
needed to describe plastic deformation of amorphous ma- 
terials near their glass temperatures; but, even there, the 
athermal analysis suggests some simplifications that may 
be useful, for example, in the choice of the transition rates 
R(s). It seems to us that the most important open ques- 
tion in the athermal STZ theory is whether the effective 
temperature x is adequate for describing all the relevant 
internal states of a deforming material, or whether other 
internal variables may be needed. Our analysis of the 
simulations by Demkowicz and Argon [T5I Ho . ITtI Il8| . 
to be described in the sequel to this paper, is aimed at 
answering this question. 
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